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Nowadays superconductors serve in numerous applications, from high-field 
magnets to ultra-sensitive detectors of radiation. Mesoscopic superconducting 
devices, i.e. those with nanoscale dimensions, are in a special position as they 
are easily driven out of equilibrium under typical operating conditions. The out- 
of-equilibrium superconductors are characterized by non-equilibrium quasipar¬ 
ticles. These extra excitations can compromise the performance of mesoscopic 
devices by introducing, e.g., leakage currents or decreased coherence times in 
quantum devices. By applying an external magnetic field, one can conveniently 
suppress or redistribute the population of excess quasiparticles. In this article we 
present an experimental demonstration and a theoretical analysis of such effec¬ 
tive control of quasiparticles, resulting in electron cooling both in the Meissner 
and vortex states of a mesoscopic superconductor. We introduce a theoretical 
model of quasiparticle dynamics which is in quantitative agreement with the 
experimental data. 

The presence of excess quasiparticles (QPs) is often characterized by an effective elec¬ 
tron temperature T which exceeds the temperature of the phonon bath Tq. The resulting 
overheating is known to be the origin of such effects as decoherence in qubit systems m. 
decrease of the quality factor of superconducting resonators BE], the excess current in 
single-electron turnstiles [B], and low efficiency of electronic cooling in normal metal (N) - 
insulator (I) - superconductor (S) junctions [3 IE]- In short, overheating is a major factor 
limiting the performance of S mesoscopic devices. More than the overall QP number Nqp, 
the critical parameter is the location of these excess QPs. For instance, for tunnel junction 
circuits, it is crucial to avoid the QPs in a superconductor nearby the junction, while the 
extra QPs located further away are of less concern. To suppress overheating in a supercon¬ 
ductor one aims at lowering the generation of extra QPs in the whole superconductor using 
proper electro-magnetic shielding of the device, and decreasing of QP density by introducing 
QP traps (see e.g. mm), by optimizing the device geometry mm, or by cooling using the 
tunnel junction to another superconductor with a larger gap HMD. The second method 
allows one to move QPs away from critical locations and relax them. QP traps have become 
an important element in designing devices for mesoscopic physics and metrology. 

The most common ones among different types of QP traps are normal metal sinks 
m, Andreev bound states in weak links |T8], special S gap engineering and 
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a) B = 0,m = 0 






FIG. 1. QP density and gap distributions in a S disc. The S gap Eg{B,x) in a S disc 
with narrow extensions (Sample A) is represented by the height of the blue volume, while the QP 
density nqp{x) is shown by red circles; m is the vorticity of the island, B is the magnetic field acting 
on the sample. The wide central part of a Sample A of size 2R is limited by vertical dashed lines, 
the narrow extensions are located on the sides, (a) Uniform zero magnetic field state; (b) Meissner 
state with reduced Eg{B,x) in the central part at small fields in a vortex free state; (c) Single 
vortex state with smaller gap reduction outside the vortex core than in (b). The 3D schematics 
depict the corresponding Eg{B,x) (in blue) and nqp{x) (in orange semitransparent). 


non-uniform superconducting states induced by an external magnetic field |H [25H27] . Here 
we focus on the magnetic field controlled trapping, a method which has a number of ad¬ 
vantages. The regions with the reduced gap in this case are of the same material as the 
rest of the device and therefore match perfectly the S parts without barriers or interface 
potentials. Besides, magnetic held gives the possibility to make tunable traps allowing, for 
example, the modulation of a resonator quality factor, needed for giant pulse formation 
(or Q-switching) in pulse lasing (see, e.g., a book [28]). The controllable use of such traps 
in various applications mentioned above assumes, certainly, understanding of their cooling 
capacities, which is necessary to optimize the designing of the particular trap conhgurations 
for different mesoscopic devices. 
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Our work aims to the solution of this ambitious and important problem focusing on both 
experimental and theoretical study of individual traps which appear in the Meissner and 
vortex states. To build a quantitative model of these traps we choose to verify it by the 
experimental measurements of the characteristics of nonequilibrium QP distributions in a 
mesoscopic S island (Al) in a single-electron transistor (SET) set-up with normal metal 
(Cu) leads. This particular device appears to provide a very convenient way to tune both 
the trap pattern applying an external magnetic held to the S island and the number of 
nonequilibrium QPs injected in the island in the Coulomb blockade conditions by operating 
it as a turnstile of single electrons [25] • The turnstile operation frequency / of the gate 
voltage modulation controls QP injection rate. This set-up allows one to probe single QP 
excitations in the superconducting dot by measuring the average turnstile current under 
pumping conditions [6] (ideally this current equals ef) and to independently control the 
vortex number in the superconductor m- The resulting trap model has perfectly proved 
its validity and efficiency in this set-up which can be used in future applications. 

Results 

Qualitative description. We illustrate the key idea of QP redistribution by Fig. [^in 
an S island with a large central part and two narrow extensions, called Sample A. In the 
absence of magnetic held acting on the sample, B = 0, the QP density Uqp is nearly uniform 
in the S island with constant gap Eg{x) = Aq, provided the heat dihusion length Lt ^ R 
is large compared to the size of the central part R (see Fig. [^a)). A small perpendicular 
magnetic held, typically few mT, which induces Meissner screening currents howing along the 
superconductor edges reduces the gap Eg[x) mostly in the wide central part of the island but 
not in the narrow extensions near the junctions [31]. Due to this non-uniform gap potential 
Eg{x), QPs illustrated by red circles are redistributed so that their density is small at the 
junctions (see Fig. [^b)). However, the total QP number is larger than that at R = 0 due 
to its exponential dependence Ngp oc minimal gap Eg^^i^ = mm^ Eg{x) 

over the island, where is the Boltzmann constant. A vortex in the island leads to further 
QP redistribution because it plays a role of a QP potential well containing a lot of QPs as 
shown in Fig. [^c). Despite its simplicity, the theoretical model that we present below yields 
a quantitative ht to the experimental data on the magnetic held and frequency dependencies 
of the pumping current, and thus to the QP distribution, rendering the turnstile an efficient 
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probe of QP dynamics and relaxation. 



FIG. 2. Layouts of the samples and results of dc measurements. Electron micrographs of 
(a) Sample A and (b) Sample B. (c) Schematic picture of the device with its electrical connections 
where the red and black lines correspond to the dc measurements, while blue and black ones to the 
pump measurements, (d) Evolution of the voltage at a fixed bias current /bias = 10 pA at the gate 
voltage CgVgje = Ug = 0.5 suppressing the Coulomb energy with the magnetic field for Sample 
A (red filled circles) and Sample B (black open squares). The vertical arrows correspond to the 
applied field values at which the vorticity m increases step by step by one from —2 to 2 as the field 
is swept from —25 mT to 25 inT. The horizontal arrow shows the direction of the field sweep. The 
solid black line is the theoretical result for this measurement. The experimental uncertainty has 
been estimated as r\j 3pV. 

DC measurements of the S gap. To probe the magnetic field induced changes in 
the gap of a S disc, we hrst measure a more basic structure, which we call Sample B (see 
Fig.j^b)). It is formed of a S disc, mimicing the central part of Sample A (Fig.j^a)), directly 
connected via tunnel junctions to normal leads at its edges. Measuring electron transport 
through the disc while applying perpendicular magnetic held H allows us to access the held 
dependence of the gap value Eg{H) at the edge of the disc and to control the vortex state. 
This way we can determine the critical helds for transitions between states with diherent 
vorticities m via simple dc transport measurements (similar approach as in Ref. [H]). We 
carried out current biased dc measurements at a gate voltage that suppresses the Coulomb 
energy (for the electrical conhguration, see red and black lines in Fig. [^c)). The experiments 
have been performed at a bath temperature of Tq ~ 60 mK (well below the S gap Aq at 
B = 0 and the Coulomb energy Ec = e‘^j{2C), where C is the total capacitance of the 
island). Note that B is the actual held seen by the sample, while the applied magnetic 
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field H differs from later due to some screening by the sample holder used for shielding the 
sample from the environment (see Supplementary Note 1 for details). 

The dc drain-source voltage V measured versus the magnetic held H, swept from —25 mT 
to 25 mT, is shown for Sample A (hlled circles) and Sample B (open squares) in Fig. |^d) at 
a hxed current of luas = 10 pA through the device. In general larger voltage corresponds to 
larger gap Eg{H) and vice versa. The sample parameters Eq, Aq, and a total normal state 
resistance across the two junctions Rt, have been extracted from IV measurements at zero 
magnetic held. For the Sample B, starting from —25 mT, the value of the voltage is small: 
the island is close to its normal state. The gap increases when decreasing the absolute value 
of the held till the maximum value reached at -1-2.5 mT with two intermediate step-like 
anomalies at H^^ —15 mT and 2 mT, corresponding to the exit of vortices, the 

hrst one from two-vortex state to one-vortex state, and the second one from one-vortex state 
to a vortex free state, respectively. Increasing H further to positive values from 2.5 mT up 
to 25 mT leads to decrease of the gap again, with two knee-like anomalies at critical held 
values 14 mT and 18 mT corresponding to the entry of the hrst and the 

second vortex, respectively. A minor distortion of the applied held (the ohset in the applied 
held 5H k, 2.5 mT corresponding to the maximal V{H) value and asymmetry of V{H) in 
the Meissner state) is caused by the sample-holder, and was corrected to theoretical curves 
only by applying the magnetization curve B{H), with B the held acting on the sample, 
measured separately (see Supplementary Note 1 for details). Note that the magnetic held 
B acting on the sample itself, is zero at the maximal V{B = 0) and corresponds to the 
symmetric V{B) = V{—B) in the Meissner state. The central part of Sample A has nearly 
the same shape and size as Sample B; thus one can expect the critical helds of these samples 
to be close to each other. The anomalies are absent in Sample A, as the gap near the tunnel 
junctions is only weakly ahected by H in the presented range. 

Theoretical analysis of DC data. For the theoretical analysis of the above exper¬ 
imental data we simplify the standard Usadel model taking into account that the size of 
the central part R of the measured samples is small compared to the characteristic length 
scale of the Green’s functions outside the vortex core regions (see Supplementary Note 2 
for details). Such approximation leads to the Usadel equation for the normal (cos6') and 
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FIG. 3. Turnstile current driven by an ac gate, (a) The pumping current I normalized by 
ef at the operating frequency / = 5 MHz, the bias voltage I4ias = 100 fjX, and at the gate offset 
rig = 0.5 versus the normalized gate amplitude Ag at applied field H = 2.4 mT {B = 0) (red dashed 
curve) and at H = —25 mT (black solid curve). The horizontal dashed lines correspond to the 
expected current quantization I = nef\ (b) the zoom-up of the first plateau at several applied field 
values between —25 mT and 2.4 mT; (c) the evolution of the current along the vertical dashed line 
in the main panel versus the field H varying from —25 mT to 25 mT; the sweep direction is shown 
by the horizontal arrow. In the further measurements we fixed the amplitude to the value shown 
by the vertical dashed line in (a). The experimental uncertainty has been estimated as ~ 10 fA, 
corresponding to 0.0125e/ at 5 MHz (not shown). 

anomalous (—isin^) Green functions 

{iE — r cos 9) sin 6* -I- A cos 6* = 0 , (1) 

with the effective depairing parameter V = ^(v^) expressed through the superfluid velocity 
= D (V(p — 2eA/hc) and averaged (..) over the sample volume (over the central part of 
Sample A) with the excluded vortex core regions. Here ip is the S order parameter phase, 
A is the vector potential determined by the magnetic held B acting on the sample, and D 
is the diffusion coefficient. The component of v^ perpendicular to the sample boundary and 
to the boundaries of vortex cores should be zero. Similarly to previous works [2H11 [25], [26] 
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the vortex cores are assumed to be normal metal cylinders of the radius Vy of the order 
coherence length ^ = ^JhD /, i.e. 0 = 0 inside the cores. The sample size 2R ~ 1 /mi is 
also smaller than the effective screening length Ae// ~ /ds ~ 2.6 /im, therefore we expect 
uniform held distribution in the island. Here A ~ 230 nm |1] is a typical bulk penetration 
depth and ds — 20 nm is the thickness of the aluminium disc. 

Solution of the Usadel equation gives us the standard expression for the hard gap Eg 
in the density of states and for the order parameter A as functions of T (see HIISHZ] or 
Supplementary Note 2). We made a ht of the held dependence of the voltage V{B) at hxed 
currents /bias using standard expressions for the current-voltage characteristic of a tunnel 
junction (see Supplementary Note 3 for details) and of the depairing parameter 

T/i\o = ai{B/Bcf-ma2B/Bc + w?a2. , ( 2 ) 

taking into account that the vector potential A in the superhuid velocity \s = D (V(p — 2eA/hc) 
is proportional to magnetic held B while the S phase distribution ip is determined by vortex 
sources. Here a* are numerical htting parameters, B^ denotes the held value of the hrst 
vortex entry and m is the total vorticity. The estimate B^ <ho/7ri?^ 10 mT based on 

^ ^ 100 nm and R ~ 0.5 /im is rather close to the value B{h2^) ~ 14.4 mT from our dc 
measurements. Here $0 = h/2e is the hux quantum. More accurate estimates of Be can be 
done numerically, e.g., within the Ginzburg-Landau approach for a concrete sample geome¬ 
try [ni[I3]. According to [7] the parameter ai determining the critical value of T/Aq for the 
hrst vortex to enter for the Usadel equation with homogeneous in a narrow strip should 
be = 0.237, while the parameters a 2 and as depend on the vortex conhguration in the 

sample. The best hts to the experimental data are obtained with ai = 0.38, 02 = 0.438, and 
as = 0.266, where we take Be = B{H^^) ~ 14.4 mT from experimental data. Parameter ai 
for a rectangular sample is expected to be a bit larger than its value in a narrow strip 

[7]. In the htting we assume that both jump-like and knee-like anomalies in the V{H) are 
associated with the change of vorticity m and verify this applying the same parameters to 
V{H) with diherent values of /bias (see Supplementary Note 3). The S gap in the narrow 
extensions of Sample A shown in Fig. |^d) is close to its zero-held value Aq up to ~ 30 mT 
with few % accuracy as the depairing parameter in this case T/Aq = {tt^wB/6 is small 
[1]. Here w ~ 130 nm is the width of the extensions. 

Pumping measurements. The pumping measurements are done on Sample A which 




FIG. 4. Excess pumping current vs magnetic field. The excess current I — ef at the first 
plateau for the driving frequencies 0.5 MHz (black), 1 MHz (red), 5 MHz (cyan), 10 MHz (violet), 
and 30 MHz (blue) with the maximum electronic temperature at B{H) = 0 of 274, 286, 320, 337, 
and 369 mK respectively. Applied magnetic field H is swept from —25 mT to 25 mT (the sweep 
direction is shown by a horizontal arrow). The vertical dashed lines show the expected values 
corresponding to the entrance of an extra vortex into the island (for H > 0) and to the removal of 
one vortex {H < 0). (inset) A close-up of the current for a field close to the expulsion of the second 
vortex. For better visibility, the data sets for different frequencies have been shifted vertically. The 
measurements were performed at the bias voltage Vbias = 100 /xV, the gate offset rig = 0.5. In 
all panels symbols (solid lines) correspond to the experimental data (theoretical model) with the 
estimated maximal QP temperature shown in numbers. The scale on the right side shows the QP 
density riqp near the junction. The experimental uncertainty is ~ 10 fA, shown as an error bar on 
both panels. 

has the highly non-uniform distribution of the gap (Fig. under magnetic field. To probe 
the magnetic held dependence of non-equilibrium QP states, we measure the current I in 
turnstile mode averaged over the period of the drive tq = 1// |29]. We apply a hxed bias 
voltage 14ias = 100 /iV and sinusoidal gate voltage CgVg/e = ng{t) = inPg + AgSm.{2'Kft) 
through the capacitor Cg with variable amplitude Ag. 

The turnstile current is expected to assume values equal to integer multiples of ef in 
the absence of non-equilibrium effects and unwanted tunneling events. The measurements 
are carried out around the gate offset point = 0.5, to maximize the expected plateau 
width, for several frequencies /. Overheating of the S island, in particular at if = 2.4 mT 
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corresponding to B{H) = 0, leads to positive deviations of this current from I = nef 
{n is an integer) by tens percents at the expected plateau positions which corresponds to 
thousands of QPs per near the junction (see Fig. |^a) at / = 5 MHz). The magnetic 
held improves QP trapping: the deviation from ef (and the corresponding QP density) at 
large enough magnetic helds decreases by an order of magnitude in the whole frequency 
range from 0.5 to 200 MHz (see Fig. and Supplementary Note 3) and approaches a few- 
percent level related to an ampliher noise, even for large gate amplitudes when pumping 
up to n = 11 electrons per cycle. The zoom up of the hrst plateau shown in Fig. [^b) 
demonstrates the magnetic held dependence of I. To separate the Meissner current from 
the vortex contribution, we present in Fig. [^c) pumping current versus the held at a hxed 
gate amplitude value indicated by the vertical dashed line in the main panel. The excess 
current 51 = I — ef increases when the held is swept from large negative values to low 
values with jumps at (see vertical dashed lines in Fig. 4). The following held increase 

to positive values leads to decreasing excess current without visible anomalies. This is related 
to the diherence in the fcth vortex entry (exit) helds Indeed, for fc = 1 at these 

helds, we have Eg{H\P) < Eg{Hp]f) which leads to the efficient redistribution of QP density 
even without any vortex (see Fig. [^b)). Despite the absence of anomalies at the vortex 
entries, it is possible to estimate the value by varying the value of the initial held: 
the discontinuous anomaly at is only visible for a held amplitude in a sweep exceeding 
hIP = 13.5 mT which is close to the value found by dc measurements in Sample B. At even 
higher values of the held, \H\ > 30 mT, the current quantization is lost again due to the 
eventual suppression of the S gap near the junctions as well (see Supplementary Figs. 4). 

Theoretical analysis of pumping data. To model theoretically the excess current as 
a function of the held B and frequency / we calculate the electronic temperature T using a 
heat balance equation 


QephiT) ~ JHbias = [ef + 5/(T)]Hbias 


(3) 


We keep in mind that T is nearly uniform and constant in time provided the heat dihusion 
length Lt is large compared to the size of the island R, and the heat relaxation time Teph, 
determined by electron-phonon coupling, is much larger than the operating period tq = 1//, 
i.e., Teph S> To allowing us to average the heat dihusion equation across the sample volume 
and over the operating period. We assume further that most of the Joule dissipation occurs 
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inside the S island and take into account the excess pumping current / averaged over the 
period from the hrst plateau on the right hand side (see Supplementary Note 5 for details). 
The value of Qeph = Q'^eph + describes the heat flow rate from the electronic subsystem 
to phonons for a non-zero depairing parameter T. We calculate its value outside 

the vortex core regions for any T combining the procedure described, e.g., in Refs. [Ml [39] 
and the solution of Eq. (§. For the experimental parameters /cbR), ^ k-^T <C Eg 

in low-temperature limit the electron-phonon heat flux 
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C(5) 


63 


k^T _j_ 


27rE4 _ 

_ ?-Te 

3ki 


(4) 


decomposes into recombination oc e and scattering terms oc e (see, e.g., [To]). Here 
E is the electron-phonon coupling constant, and V (V^,) is the volume of the island (the vortex 
core regions). In the vortex cores in the same limit of negligible phonon temperature Tq T 
the electron-phonon heat flow is modelled by the standard normal metal expression = 

SVuT® with the volume of m vortex cores assumed to be Vv = mrlds- The recombination 
term in Eq. ( [I7j ) becomes dominant at k-^T > O.lEg. Beyond the low temperature limit 
we use a numerically calculated expression for Q^XiT) instead of (17) (see Supplementary 
Note 5 for calculation details). 

Eventually we obtain the magnetic held and frequency dependence of the measured excess 
current SI{T) as 


6I{T) = C 


y'2'KAokBT 

bRt 


exp[-Ao/kBT] 


(5) 


Note that the QP density near the junction is proportional to the excess current Ugp = 
D{Ep)eRT5I{T)/C and can be extracted from SI(T) using the normal state density of 
states in the superconductor D{Ep) (see the scale on the right side of Fig. [^showing the 
QP density riqp). Here C ~ 1 is a numerical coefficient determined by the wave-form and 
the amplitude Ag of the gate drive, in particular the duration for one junction to be open 
for tunneling in each cycle. A detailed derivation is given in Supplementary Note 6. Note 
that in Eq. ([^ we neglected contributions of higher order processes in R^ like Andreev 
tunneling due to the small transparency of the junctions (see experimental results in nzi 
and estimates in Supplementary Note 6). By solving (|^ with the substituted expressions (|^ 
we hnd the solution for T and SI{T) (solid lines in Fig. |^. We used the constant C = 1 
for a hxed drive amplitude Ag. The main uncertainty in the htting procedure originates 
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from the parameter EV. The volume of the S sample can be estimated based on the electron 
micrograph (see Fig. [^b)) as V ~ 3 ■ 10“^° m^, but usually this value is overestimated due 
to additional uncontrolled oxidation of Ah On the other hand the typical range of the mea¬ 
sured values of the electron-phonon relaxation constant S in the bulk aluminium [HI [15], Il2] 
is within 2 ■ 10® to 5 ■ 10® W m“®. Our htting gives results agreeing reasonably well 
with the experimental data within the range of SV from 4 ■ 10“^^ to 9 ■ 10“^^ W K“®. In 
Fig. I^we present a £t for a certain middle value EV = 6 ■ 10“^^ W which is in the best 
agreement with the experiment at moderate frequencies. Assuming V ~ 3 ■ 10“^° m® we get 
E = 2 ■ 10® W m“® which is towards the low end due to the overestimated V but within 
the range given above. We have extracted the optimal value of the vortex core radius within 
the range = 2.5 — 2.7^ both from the dc measurements (see Supplementary Note 3) and 
from the pumping data, which is in perfect agreement with the previous theoretical results 

[2113]. 

Discussion 

According to the theoretical model, Eq. ([^, the maximal electronic temperature at Ag = 
1.1 and / = 30 MHz is T ~ 370 mK. It corresponds to a number of non-equilibrium QPs 
Ngp = UgpV ~ 250 in the uniform state (see Fig.j^a)). In the held increasing from B{H) = 0 
the Meissner supercurrents sufficiently improve the electron-phonon relaxation by reducing 
the gap Eg{r) in the central part of the island even before the hrst vortex enters the island. 
This leads to the at least 10—20 times reduction of the QP density near the junction when the 
excess current approaches the ampliher noise level. The vortex contribution is clearly seen 
in the decreasing held regime due to the hysteresis caused by vortices. Indeed, the vortices 
that entered the island at a certain value of the held stay there till smaller helds (where the 
ehect of Meissner current is smaller) and improve the relaxation of hot QPs most ehectively. 
Such hysteresis allows us to see the vortex contribution alone (see the larger step in Fig. at 
H ~ —2 niT) and the improvement of relaxation in the two-vortex state with respect to the 
one-vortex state (the smaller step at 77 ~ —13 mT). We estimate the recombination rate 
in the vortex state Tree — f /Ngp^vort as the injection QP rate / divided by the QP number 
Ngp^vort — 2D{Ep)Vvk-BT\\i2 in the vortex core volume Vy (see Supplementary Note 5 for 
details). At / = 30 MHz it gives Ngp^yort ~ 100, F^ec — 0.3 MHz of the recombination rate, 
i.e., 20 times higher than F°g^ = 16 kHz estimated in |T5| at 7? = 0. 
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In conclusion, we demonstrate effective control of the number of excess quasiparticles 
and their spatial distribution in a mesoscopic superconducting disc by applying a small 
magnetic held on it. We hnd that both the Meissner supercurrents and vortices entering 
the disc one by one each give important observable contributions to the trapping of non¬ 
equilibrium quasiparticles. We demonstrate that a single-vortex contribution is sufficient 
to keep the superconducting disc near equilibrium up to 30 MHz injection frequency with 
Uqp ~ 400 /im“^ quasiparticle density near the junction and recombination rate of order of 
Tj-ec — 0.3 MHz. Our dc and pumping measurements conhrm the assumption Eli that a 
vortex can be considered as a normal metal cylinder with the effective radius = 2.5 — 2.7^ 
both in charge and heat transport problems. Our theoretical analysis of the quasiparticle 
trapping has proven its validity and efficiency in the set-up being in quantitative agreement 
with the experimental data. 


Methods 

Device fabrication. The hybrid devices with aluminium as the superconductor, copper 
as the normal metal, and aluminium oxide as the tunnel barrier in between, have been fabri¬ 
cated by standard electron-beam lithography and two-angle shadow evaporation technique. 
The aluminium island is ds = 20 nm thick and it is oxidized with O 2 for 2 min at 2 mbar. 
The copper leads, 25 nm thick, are placed on the oxidized A1 forming tunnel junctions. 

Sample geometries and parameters. Two different island geometries have been 
employed in the measurements: Sample B has a nearly square-shaped island, as shown 
Fig. [^b), and Sample A with the same central part as geometry B has two additional 
long narrow aluminium extensions from each side toward the junctions (Fig. [^a)). The 
diagonals of the island are 2R ~1 pm both in A and B, and the narrow extensions of the 
island in A are 2 pm long and w ~ 0.13 pm wide. The sum of the tunnel resistances 
of the two junctions is Rt — 577 kH for Sample A and Rt — 714 kH for Sample B. We 
measured the IV characteristics of SETs at various values of the DC gate voltage at the base 
temperature to determine the zero held S gap value Aq — 190 (207) peV and the charging 
energy Eq ^ 173 (133) peV for Sample A (B). 

Reproducibility and noise. All the results presented here are reproducible between 
different runs and between samples of similar geometry, in particular, as concerns the values 
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of the critical fields of vortex entry (exit). The results depend only on whether the absolute 
held value increases or decreases, provided by the hysteresis in vortex entry/exit events, but 
they do not depend on the sign of the held as such. The samples are cooled down through 
the superconducting transition with a zero-held cooled magnet. The uncertainties of current 
and voltage are estimated to be 10 fA and 3 /rV respectively. They are taken as the noise 
from the amplihers. 


* Present address: LPMMC, CNRS/Foundation Nanosciences under the aegis of Joseph Fourier 
University Foundation, BP 166, 38042 Grenoble, France 
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Supplementary Materials 


SUPPLEMENTARY NOTE 1. HYSTERESIS IN THE MEASUREMENTS UN¬ 
DER FIELD 

As pointed out in the main text, a remanent field of 6H ^ 2.5 mT and the asymmetry 
in the Meissner state are present at the sample location at zero applied field H = 0 due to 
the presence of superconducting parts in the sample-holder. The dc measurement presented 
Fig. 2(d) of the main text for Sample B has been reproduced (on the same sample) in 
a sample-holder that does not have the superconducting shield for which B = H. These 
measurements are shown in Supplementary Figure]^ for several bias current values /bias = 1 
(blue triangles), 10 (red circles), and 100 pA (black squares), together with the theoretical 
model for the experimental data (solid lines of corresponding colors). The dc voltage V 
vs field B in this measurement is symmetric with respect to the zero applied field value 
(except for the vortex hysteresis intrinsic for the sample). The theoretical model presented 
in Supplementary Note 3 reproduces perfectly the experimental points. 



FIG. 5. DC measurements without field distortion. Evolution of the voltage at a fixed bias 
currents /bias = 1 (green triangles), 10 (red circles), and 100 pA (black squares) at the gate voltage 
CgVgle = ng = 0.5 suppressing the Coulomb energy with the magnetic field for Sample B in a 
sample-holder that does not have the superconducting shield together with the theoretical model 
(solid lines of corresponding colors). The field B is swept from —30 mT to 30 mT. 

In order to fit the theoretical model to our measurements versus field H performed in the 
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sample-holder with the superconducting shield (i.e. with deformation of the held prohle), a 
correction has to be done, as the held H applied through the coil dihers from the ehective 
“acting” held B seen by the sample. The held prohle has been measured in the shielded 
sample-holder with a Hall sensor at 4.2 K and at 0.2 K. The hrst measurement has been 
done in the normal state when there is no magnetic shielding as a reference point (not 
shown). The second measurement of the ehective held B versus the applied coil magnetic 
held H swept from —30 mT to 30 mT (at 0.2 K) is shown in Supplementary Figure |^as black 
dashed lines. The arrows point out the direction of the sweep. For \H\ > 20 mT, the sample- 
holder is fully normal and the ehective held equals the applied one. For \H\ < 20 mT, a 
nonlinear superconducting response from the sample holder is present leading to hysteresis. 
The red line in Supplementary Figure is the correction found by comparing directly the 
dc measurements in both sample-holders {V{H) in Fig. 2(d) of the main text and V{B) in 
Supplementary Figure]^. The two methods yield very similar results. Eventually theoretical 
curves given as functions of B by the model are presented as functions of H using red curve 
B{H) in Supplementary Figure]^ 



FIG. 6. Magnetization cnrves. Field correction found by using a Hall sensor (black dashed 
lines) and by comparing the dc measurements of V(H) in both sample holders (red solid line). The 
largest deviation from linearity lies at low field. The two methods yield a very similar correction. 
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SUPPLEMENTARY NOTE 2. HOMOGENEOUS APPROXIMATION IN A 


MESOSCOPIC SAMPLE 


The samples studied experimentally are in the dirty regime, namely i ^ where ^ = 
^hD/Ao is the coherence length, D is the diffusion coefficient, i is the elastic mean free 
path, and Aq is the superconducting gap. The quasiparticle (QP) spectral characteristics in 
this case can be found from the Usadel equations (see, e.g., m) 


hD„r,^, , / „ hr. . 

—V 0(r) + UE - —V, cos0(r 


sin6*(r) + Acos6'(r) 


0 , div (sin^ 6 '(r)v 5 ) = 0 , (6) 


where = —g^* = cos0(r) and = —f"^* = —fsin6'(r) are the normal and anomalous 
Green’s functions (superscripts ’R’ and ’A’ stand for ’retarded’ and ’advanced’). 

Considering the experimental situation of a mesoscopic superconducting sample with the 
characteristic size R = 0.5 /im and the coherence length ^ ~ 100 — 200 nm (depending on the 
diffusion coefficient), we have to verify if we can neglect the gradient terms in Sup. Eq. (|^. 
For subgap energies the characteristic length scale of the function 6*(r) can be estimated 
as follows: — E^/A^. It is natural to assume the ^(r) function inhomogeneity to be 

small provided fy -^/l — E^/A^ > R. This condition gives us the energy interval 1 — E‘^/A'^ < 
i^/R^ ~ 0.1, sufficient for the calculations of the electron-phonon heat flow Qeph and the 
thermal excitation leakage current 51 for temperatures much lower than the superconducting 
gap. Indeed, the main contribution to Qeph and to 51 is given by |E/A — 1| ~ k^T/A ~ 
0.03 < ^ Q i 

Certainly the above assumption is strictly valid only for the Meissner state: inside the 
vortex core the gap and the anomalous Green function turn to zero at the scale of the 
effective core radius which is of order of the coherence length ^ [2] . 

To avoid numerical solution of the Usadel equation we adopt in the main text the following 
approximate procedure. In the presence of vortices we assume that both the order parameter 
and 6{r) function vanish inside the vortex cores while outside the core regions we assume 
the 6*(r) function to vary slowly and introduce, thus, its average 9 over the region outside 
the vortex cores (omitting the spatial dependence in the notation). The deviations from the 
averaged order parameter A beyond the cores also become small in this limit. Integrating 
now the above Usadel equation over the region outside the vortex cores we obtain Eq. (1) 
from the main text with the effective depairing parameter expressed through the superfluid 
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velocity as 


r = ^(v^) = ^((Vv7-2eA/nc)^) . (7) 

Here the brackets (..) denote an average over the sample volume (over the central part of the 
sample A) with the excluded vortex core regions, (p is the superconducting order parameter 
phase and A is the vector potential determined by the external magnetic held B applied to 
the sample. The second Usadel equation in our approximation reduces to divv^ = 0 and 
leads to the vanishing components of v* perpendicular to the sample boundary and to the 
boundaries of vortex cores. Here and further on we neglect the changes in the magnetic 
held B due to the screening currents howing in the sample due to the smallness of the 
characteristic sample size R as compared to the ehective screening length Ae// = X'^/ds- For 
our samples A ~ 230 nm |1] and ds = 20 nm, therefore Ag// — 2.6 /rm. 

Solution of the averaged Usadel equation gives us the expression for the hard gap Eg in 
the density of states and for the order parameter A as functions of F as PIHHT] 

Eg = A(1 - 72 / 3 ) 3/2 ^ ^ _ Aoe-’^^/^ , 7 = F/A . ( 8 ) 

In the main text we focus on the case 7 < 1 (F < implying that the gap Eg > 0 

is non-zero. 

SUPPLEMENTARY NOTE 3. DC FITTING 

Using the solution of the averaged Usadel equation, Eq. (1) from the main text, one 
can ht the IV characteristics shown in Supplementary Figure Indeed, we consider a 
hybrid single electron transistor (SET), namely, a mesoscopic superconducting island tunnel 
coupled to the normal metal leads (NISIN). We apply a hxed bias current /bias through the 
normal leads and the constant gate voltage Ug = CgVg/e = 0.5 to the gate electrode coupled 
to the island through the capacitor Cg (see black and red lines in Fig. 2(c) of the main text) 
and measure the difference V of voltages Vl^r = ±U/2 applied to the leads as a function of 
the magnetic held B seen by the sample. 

In stationary state the current /bias flowing from one lead to another is equal in any cross 
section and it can be calculated in any of two junctions (for example, in the left one) 

k 
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as a sum over the island charge state k of the sequential tunneling rates to 

(from) the island through the left junction. This sum is weighted with the probability pk of 
system being in this charge state, which is calculated using the standard rate equation for 
the balance of the probability fluxes gHin] 

k±l^kPk±l — ^k^k±lPk] , '^^Pk = 1 ; ( 10 ) 

in the stationary case dpk/dt = 0 with the tunneling rates Ffc-s-feii = ^ 1 , 1=1 R^''k^k±i ^^0 
^Uk±i = Wt] given by 

r[(7] = ^/ ns{E)ME)ll-fT„{E+U)]dE . (11) 

Here = ^2Ec{k—ng±l/2)EfeVi are the energies gained by the electron tunneling to/from 
the island (being in the charge state k) through ith junction, Rt/2 is the tunnel resistance 
of each junction. Here we focus on the magnetic held effects in the sample B (see Fig. 2(b) 
in the main text) and neglect all the overheating effects assuming the equilibrium Fermi 
distribution of electrons over energy fxiE) = + 1 ]“^ with the electron temperature 

T equal to the phonon bath temperature Tq. The density of states (DOS) ns{E) = Re[cos 6 '] 
normalized to its normal state value D{Ef) in the superconducting (S) island near the 
junction is obtained from the solution of averaged Usadel equation, Eq. (1) from the main 
text, with the depairing parameter F/Aq = ai {B/Bc)‘^ — ma 2 B/Bc + Eq. (2) in the 

main text, having three positive numerical htting parameters a;. 

In the Sample B the tails of the wave functions localized in the vortex core(s) give a 
substantial contribution to the DOS and to IV curves for m 7 ^ 0 at small /bias, but they 
are not included into the averaged model. To model this contribution we replace the DOS 
ns{E) by ns{E){l — ]^y adding the phenomenological normal metal DOS 

with the exponentially suppressed prefactor determined by the vortex distance from 

the junction R ~ 0.5 /im and by the exponential decay of the wave function localized in the 
vortex core of the effective radius r^. This vortex contribution leads to reduced V{B) at 
small /bias in the mixed state m 7 ^ 0 and to the suppression of the jumps at the vortex entry 
helds (see red and green curves in Supplementary Figure]^. 

By htting V{B) at /bias = 100 pA which is not affected by the vortex tail contributions 
one can extract the following values of htting parameters ai = 0.38, 02 = 0.438, and 
Q !3 = 0.266 mentioned in the main text. Following m we attribute to all jumps in this plot 
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with the change of the number of vortices in the sample and use the point of the hrst jump 
at -B > 0 as the held of the hrst vortex entry = 14.4 mT. In this setup we don’t see 
any transitions between vortex conhgurations with the constant vorticity like the transition 
to a giant vortex state (see, e.g., [El [13]). Using these parameters one can ht V{B) quite 
well at all bias current values with R/r^ ~ 1.7. The optimal value of the vortex core radius 
= 2.5 —2.7^ extracted from dc measurements in the Sample B is in perfect agreement with 
the previous theoretical works |2], |3]. In subgap regime /bias = 1 and 10 pA the jump-like 
anomalies in V{B) become knee-like, but because of the above-mentioned reasons we still 
associate each of them with the vortex entry or exit. 


SUPPLEMENTARY NOTE 4. ELECTRONIC PUMPING 

The electronic pumping of the Sample B at / = 5 MHz, when Vbias ~ 120 /iV and 
Ug ~ 0.5, is shown in Supplementary Figure [^a) with the held H swept from —10 mT 
to 2 mT. Contrary to what is observed in the Sample A, the increase of the magnetic 
held increases the deviation from the current quantization I = ef, due to the ehect of the 
screening current on the superconducting gap. This observation is in agreement with the 
theoretical model with the increasing number of QPs in S island with the held. 

The electronic pumping of the Sample A at frequency / = 200 MHz, when Vbias ~ 250 /iV 
and ng ~ 0.5, is shown in Supplementary Figurej^b) for two held values H = 0 and —20 mT. 
Similarly to the lower frequency range the increase of the magnetic held reduces the QP 
contribution to the excess current. 

The evolution of the pumping current in the Sample A with the held is shown in Supple¬ 
mentary Figure]^ In panel (a) we show the pumping current versus held with the diherent 
initial held values. For a small value of the initial held H = —12 mT, the island is in the 
Meissner state and a continuous variation of current is observed. The anomaly at —2 mT 
appears only if the initial held is large enough to turn island into the mixed state (see blue 
and red curves). Similarly, the anomaly at ~ —15 mT (see inset of Fig. 4 of main text) 
appears only if the initial held exceeds 20 mT. Panel (b) shows that the current deviates 
signihcantly from I = ef at helds larger than ~ 30 mT due to the reduction of the S gap 
near the junctions. 
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FIG. 7. Pumping of Sample B and sample A at 200 MHz. (a) Electronic pumping of the 
Sample B at / = 5 MHz at fixed bias voltage Ihias = 120 /iV under field H from —10 mT to 2 inT 
(from top to bottom). With the field, the deviation from I = ef increases rather than decreases 
in contrast to the behavior of Sample A. (b) Electronic pumping of the Sample A at / = 200 MHz 
at fixed bias voltage 14ias = 250 /xV in field H = 0 (red dashed line) and H = —20 mT (solid black 
line). 

SUPPLEMENTARY NOTE 5. HEAT BALANCE EQUATION 

In this section we describe the theoretical model of the relaxation of QPs in applied mag¬ 
netic held by using the example of NISIN SET in the turnstile regime. By applying the 
constant bias voltages Vl,r = ±14ias/2 to the normal leads and the periodic gate voltage 
ng{t) = CgVg{t)/e = n^g + AgSini^nft) with a certain offset n°, frequency /, and the ampli¬ 
tude Ag to the gate electrode one can push electrons to tunnel through the system producing 
a time-dependent current I{t). This transport current I{f) howing from one lead to another 
drives the NISIN turnstile out of the equilibrium by injecting nonequilibrium QPs into the S 
island. The power Q%is injected to the island increases with the frequency / and we model 
this increase in mean density of QPs in the superconductor by raising its electron tempera¬ 
ture T relatively to the phonon bath temperature Tq. Note that the quasiequilibrium Fermi 
distribution of electrons over energy /t(E) = jg provided by the smallness of 

the inelastic electron-electron scattering time Tee comparing to the operating time tq = 1 // 
and the effective charging time e//.|14j 

Due to the large electron-phonon relaxation length Lt ^ R compared to the island size 
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FIG. 8. Extra pumping measurements of Sample A. (a) Evolution of the electronic pumping 
with the magnetic field measured in the range ±12 mT (black squares), ±15 mT (blue open dries) 
and ±25 mT (red crosses). At the initial field value the island is in a vortex free, in a single-vortex, 
or in a two-vortex state, respectively, (b) Pumping current in a wide field range, from —60 mT 
to 60 mT. The pumping has been measured with the gate offset rig = 0.5 and the bias voltage 
Gbias — 100 /uV at / = 5 MHz. 


R we consider the heat balance equation |1] for the whole island 


Q NIS Qeph ) 


( 12 ) 


where Qeph is the electron-phonon heat flow averaged over the island volume V. The power 
injected into the island from the junction can be written as follows 


QnIS = Ep* + Qs(K,) 

k^i 

as a sum of the power dissipated in the S island in a single tunneling event 

1 


(13) 


Qs{U) = 


e^Ri 


Esns{E)fT{E) [1 - fTg{E ± U)] dE , 


(14) 


over the junction i = L,R through which electron tunnels and over the island charge state 
k. This sum is weighted with the probability pk of system being in this charge state, which 
is calculated using the standard rate equation for the balance of the probability fluxes Sup. 


Eq. (10). Here Eg = E. Note that here and further on we neglect the difference between 


the electronic temperature T^r in the normal metal of volume Vat and the bath temperature 
To, because of sufficient electron-phonon relaxation there Q^ph = — Tq) (Ejv is 
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the electron-phonon material constant) and the rather small power Q^js injected into the 


normal leads which can be obtained from Sup. Eq. (13) by replacing Es by = eVi — E 


in Sup. Eq. (14). 


Due to the Joule’s law the sum of the powers injected into the island Qffjs into the 
normal leads Q%js should be equal to IV, where / is the current in Sup. Eq. 0 averaged 
over the period tq = 1// of the gate voltage operation. Usually if the superconductor is 
not overheated much k^T Eg one can assume that Qf^js 3> Qjy/s subgap regime and 
Qffjg ~ Qnis biases above the gap [Hj. As a result within constant factor of order of 
unity the heat balance equation can be rewritten as follows 


Qeph — IV . 


(16) 


Within the optimal conditions of the proper turnstile shielding and optimized device geom¬ 
etry the averaged current I is close to its ideal value ef and the deviation 51 = I — ef is 
mainly governed by nonequilibrium QP density in the S island near the junction 


n,p(r) = 2D{Ef) / ns(E)fT{E)dE a D{Ef) . 


(16) 


The estimates for contributions in higher orders in small parameter h/ c^Rt are given in the 
following Note. The latter expression is written for the sample A where the S gap near the 
junction is close to Aq. At 5 = 0 the QP number Ngp in the S island equals Ngp = UqpV the 
product of Ugp to the volume of the island V. The QP number Ngp^^ort — 2D{Ep)VykBT ln2 


in the vortex core can be obtained from Sup. Eq. (16) by substitution of the normal state 


DOS ns{E) = 1 into the integral and by multiplying it by the vortex core volume Vy = 
Tyds ~ 7^‘^ds. 

The electron-phonon heat flux under magnetic held is similar to the expression given by 
Eq. (3) in [T5] 

SV 


Qeph 


[nrie) - riTgie)] de / Me^e+e [f{E) - f{E + e)] dE 


(17) 


24C(5)fcB Jo j-oo 

with the term ns{E)ns{E -f e)[l — Ag/E(E -|- e)] substituted by ME,E+e of the form [T6] 

8Me,e' = 2{gE' — gE')i9E ~ 9 e) ~ (Je' ~ ~ fl; ) ~ Ue ~ fE)ifE^ ~ Je') ■ (^ 8 ) 

Here S is the electron-phonon material constant, and C(s) is the Riemann zeta function. The 
retarded (advanced) normal and anomalous f^V) Green’s functions are determined 
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by the solution of Eq. (1) from the main text, i.e., Me,e' = ns{E)ns{E') — b{E)b{E'), with 
ns{E) = Re[cos6'] and b{E) = Im[sin6']. 


In the low temperature limit Tq,T Eg/k-Q the main contribution to Sup. Eq. (17) 
arises from the energies 0 < \E\ — Eg < k^T close to the hard gap value EEg, which can be 
calculated using the following expansion of ns{E) and b{E) over the small positive parameter 
5E = \E\ 


Eg^Eg 


ns{Ef = Q{5E) 


2(5EA2/3 


b{E) 


if 


(19) 


’ ns(E)s\gn(E) 

Here Q{x) is the Heaviside theta-function. 

Substituting Sup. Eqs. (18, 19) into Sup. Eq. (0 and taking into account only the 
leading terms in the small parameter k-^T/Eg we obtain 


Qeph 


sv 


64^, - 
—T^e 


+ 


2TTE^g 


-Te" 


( 20 ) 


C(5) i63' 

for rather large electronic temperatures To ^ T -C Eg/k-Q. Note that the 

recombination term (cx dominates at k-^T > Q.lEg and should be taken into 

account. In the opposite case Tq <^T /k^ 

SVT^ r 128 


Qeph 


n 'XEl '2Eg 

— k^T e .eT 


( 21 ) 


9C(5)r2/3£:j''‘‘ 

Here the recombination term is of order of the scattering term (cx at ksT rs-/ 0.3Eg. 

In both cases as the temperature becomes of the order of the gap one have to use full 
numerical expression given by Sup. Eq. 


SUPPLEMENTARY NOTE 6. EXCESS CURRENT AS A FUNCTION OF ELEC¬ 
TRONIC TEMPERATURE 


To calculate our main observable, the leakage current SI = I—ef in the NISIN turnstile we 


use the simplihed version of the master equation given in Sup. Eq. (10) for low temperatures 
taking into account only two charge states k = 0 and k = 1 


^ = To^iPo - hi^oPi , Po = 1 - Pi • 


( 22 ) 


with the tunneling rates given in Sup. Eq. (0 in the subgap regime < Eg given by 

r[U] ^ (23) 













Here = 2EcAgSin{2TTft) — eVi and Tt = ^27rk-BTEg/e^RT- This expression 

contains the exponentially growing part with U which determines the dominant tnnneling 
rate with maximal U for each time instant. 

We consider the offset = 0.5 for simplicity and use the symmetry of the drive ng{TQ — 
t) = 1 — ng{t) focusing on the first half of the period with Ug increasing from 0.5 — Ag to 
0.5 + Ag. We assume that before the time instant ti the island is discharged k = 0 due to 
the domination of rate rf_^Q among the others and the charging process is started ai t = ti. 
The probability Poit) to stay in the state fc = 0 is decreasing with time t > ti as 

pt 


Po(t) = exp 


- / To^,{t')dt' 


(24) 


' ti 


For typical frequencies the charging process occurs not far from Ug = 0.5, therefore further 
we linearize the drive ng{t) ^ 0.5 + 27rAg{ft — 1/4). As the island has been charged Po{t*) = 
e < 1 (let’s take e = 1/2 for definiteness) the leakage current starts to flow. The number of 
excess electrons Ni through the island can be written as the integral of the largest subleading 
rate rf_^Q(f) governing the leakage current over the time interval t* < t < t 2 before this rate 
becomes the dominant one 


f*t2 


N, ~ 


1-s-O 


{t)dt 


(25) 


The leakage current can be calculated as follows 61 ~ 2efNi, where ’2’ accounts for the 
leakage during the second half of the period due to the symmetry k -H- 1 — k and L -H- P. 


By substituting Sup. Eq. (23) in Sup. Eqs. (24, 25) and calculating the time instants 
ti = Tq — t 2 and t* one can come to the result 


SI = I-ef eTj 


2Eg - |e|l/ - ksToaiT) 


(26) 


2'KEcAg 

where a = In 27 rE^. 4 ^Jin 2 %(^i) = %(^ 2 ) = 0.5 at low enough electronic temperature 

T < Toll - (eVI2E,)]-^ and a(T) = |e|V72 - ||(1 - r„/T) + In for n„(U) = 

1 — ng{t 2 ) < 0.5 in the opposite case T > To[l — {eV/2Eg)]~^. In this derivation we consider 
the operating frequency / to be small compared to the charging rate ~ (f* — fi)“^ to avoid 


missing 


events. We neglect the relative corrections of order of e and 


rf_,,o/rf-s.o when the first term in Sup. Eq. (23) dominates for all rates). We 

don’t take into account the factor 1/2 in W during the time when S> rfy.g = rf_,^Q ~ 
when the discharging occurs with the equal probability Pl,r = 1/2 to the left 
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and to the right contact. We can do it, becanse dnring the integration of Snp. Eq. (25) 


can go beyond the snbgap range < —Eg snppressing the second term in Snp. Eq. (23) 
for exponentially ~ and keeping the rate rf_^Q to be the dominant one 

in the leakage current. 

To avoid all these unimportant details we consider a certain A^-dependent numerical 


prefactor C ~ 1 instead of the square brackets in Sup. Eq. (26) and come to Eq. (5) of the 
main text by using the assumption that the S gap near the junction (in the sample A) is 
close to Aq. 


Comparing Sup. Eq. (16) and Eq. (5) in the main text one can write down the following 


relation between the leakage current 51 = I — ef and the QP density n™ near the junctions 


Hqp = D{EF)eRT5I/C 


(27) 


used in Fig. 4 of the main text to show the QP density scale. 

Note that we neglect also the contributions of higher orders in the small parameter 
h/^e^Rr) like Andreev tunneling (see, e.g., [El HI]) due to rather large tunnel resistance of 
the sample contacts. Indeed, from the experimental side the attribute feature of Andreev 
tunneling is the additional peak in the beginning of each current plateau I = nef ini 
which is not observed in all pumping measurements of this paper. From the theoretical side 
one can estimate the relative contribution 51 ar/{ ef) of Andreev tunneling to the current 
as the ratio Pyi_R/r[f/] of dc rates of sequential T[U] = UKc^Rt) and Andreev tunneling 
PyiR ~ TThU / {AN Rff in the above-gap regime. Here N = A/Ach is the number of channels 
in the tunnel junction, A ~ 6 ■ 10^ nm^ is the area of the junction and Ach is the are of a 
single channel. Theoretical estimates given in [18] lead to A^h ~ 2 nm^, while experimental 
observation [E| gives Ach 30 nm^. The upper bound estimate with A^h ~ 30 nm^ and 
Rt = 577 kH for the sample A gives N ~ 200 and 5lAR/{ef) ~ 7ih/{ANe^RT) — 3 ■ 10“® 
which can be neglected comparing to the QP contribution. 


* Present address: LPMMC, CNRS/Foundation Nanosciences under the aegis of Joseph Fourier 
University Foundation, BP 166, 38042 Grenoble, France 
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